Hysteresis loops and adiabatic Landau-Zener-Stiickelberg 
transitions in the magnetic molecule {V 6 } 
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We have observed hysteresis loops and abrupt magnetization steps in the magnetic molecule {Ve}, 
where each molecule comprises a pair of identical spin triangles, in the temperature range 1-5 K for 
external magnetic fields B with sweep rates of several Tesla/ms executing a variety of closed cycles. 
The hysteresis loops are accurately reproduced using a generalization of the Bloch equation based 
on direct one-phonon transitions between the instantaneous Zeeman-split levels of the ground state 
(an S = 1/2 doublet) of each spin triangle. The magnetization steps occur for B as and they are 
explained in terms of adiabatic Landau-Zener-Stiickelberg transitions between the lowest magnetic 
energy levels as modified by inter-triangle anisotropic exchange of order 0.4 K. 

PACS numbers: 75.50.Xx, 75.45.+j, 71.70.-d 



Magnetic molecules provide a very convenient plat- 
form for exploring fundamental issues in nanomagnetism. 
Heisenberg exchange between the magnetic-ion spins em- 
bedded in each molecule gives rise to a discrete spectrum 
of magnetic energy levels. Moreover, the magnetic inter- 
action (dipole-dipole) between molecules is generally so 
small as compared to intra-molecular exchange interac- 
tions that a crystal sample may be regarded as a macro- 
scopic assembly of independent identical quantum nano- 
magnets. One significant goal is to understand the inter- 
actions of the magnetic molecules with the environment 
("heat bath"), for example via phonons. In particular, 
it is essential to understand the nature of the thermal 
relaxation mechanism, the controlling factors responsi- 
ble for irreversible and dissipative phenomena, and the 
detailed route to thermal equilibrium of these nano-size 
quantum spin systems. The simpler the spin system the 
greater the prospects for achieving a deep understanding 
of the underlying issues, and this opportunity is provided 
by the magnetic molecule {Ve}0. Each {V$} includes 
a pair of triangles of exchange-coupled vanadyl (V0 2+ , 
spin 1/2) ions. As shown below, at low temperatures the 
instantaneous magnetization, M(t), of this spin system, 
in response to pulsed magnetic fields, B(t), with sweep 
rates of several Tesla/ms, exhibits pronounced hystere- 
sis loops as well as abrupt magnetization steps that are 
due to Landau-Zener-Stiickelberg (LZS) transitions^ 
between lowest energy levels. By explaining the de- 
tails of the dynamical magnetization one establishes both 
the low-temperature relaxation mechanism for the indi- 
vidual magnetic molecules as well as microscopic infor- 
mation concerning the lowest energy levels, not readily 
accessible. Indeed, our analysis suggests the existence 
in this magnetic molecule of an effective inter-triangle 
anisotropic exchange of order 0.4 K; otherwise Kramers' 
theorem^, [j| would forbid the occurrence of LZS transi- 
tions. 



There are several important differences between the 
present work and previous studies of M(t) in magnetic 
molecules in time-dependent magnetic fields. From our 
observation of hysteresis effects in {Vq} we conclude that 
the thermal relaxation time r in this molecule is of or- 
der 0.1 ms. This is many orders of magnitude shorter 
than those reported for "single-molecule magnets" such 
as {Mni 2 }0 and {Fes}0 where a large anisotropy en- 
ergy barrier is responsible for relaxation times of order 
10 3 — 10 5 sec. Also, we assume that the phonon bottle- 
neck effect which typically occurs at low temperatures 
(e.g. T < 200 mK for {Vi 5 As 6 }[l]) does not arise: For 
the temperatures of our experiment (T > 1.5 K) the num- 
ber of available resonant phonons per molecule is large so 
that they equilibrate independently from the spins (typ- 
ical times T p h < 10~ 6 s, much smaller than both the 
experimental time scale T exp ~ 1 ms and the relaxation 
times r of the spins) . Moreover, due to the high sweep 
rate of B(t) in our measurements, LZS transitions are 
consequential only in the immediate vicinity of B = 0. 
Away from B = 0, we use a generalization of the stan- 
dard Bloch equation for M(t), where the relaxation rate 
depends on the instantaneous B(t). The excellent agree- 
ment obtained between theory and experiment allows us 
to identify the dominant mechanism for thermal relax- 
ation in terms of direct one-phonon processes. To our 
knowledge, this is the first time that quantitative agree- 
ment between theory and experiment has been achieved 
for hysteresis loops in magnetic molecules. 

We first summarise the most important known features 

of {V 6 }Q. The magnetic molecule [H 4 vJ V 8 (P0 4 )4 
{(OCH 2 ) 3 CCH 2 OH} 2 ] 6 -, abbreviated as {V 6 }, and iso- 
lated as (CN 3 H 6 )4Na 2 {V6} • 14 H 2 0, may be pictured 
(see Fig. 1) in terms of two identical triangular units 
per molecule, each unit consisting of three spins s = 1/2 
(V0 2+ ions) interacting via isotropic antiferromagnetic 
(AFM) exchange. Two of the 2-spin exchange constants 
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pendent S — 1/2 spin triangles. 
Energy (a.u.) Time-resolved magnetization measurements were per- 




formed on a powdered sample for half-cycle and full-cycle 
sweeps by a standard induction method using compen- 
sated pickup coils and a nondestructive long pulse mag- 
net installed at ISSP. Utilizing fast digitizers, the induc- 
tive method provides data for dM/dt and dB/dt which 
are subsequently integrated to give results for M versus 
B. The pulsed fields have a nearly sinusoidal shape as 
a function of time, with a half-period about 21 ms (0- 
maximum-0). The sample of 36.9 mg was packed in a 
thin-walled cylindrical teflon capsule (inner diameter 3.0 



-6 -4 -2 2 4 6 mm) and then directly immersed in a liquid Helium bath. 
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FIG. 1: Left: Structure of the two spin triangles in the 
{Ve} anion (bright grey spheres). Phosphate exchange paths 
(green), and other ligands not shown, mediate strong intra- 
triangle exchange (J a , blue) and weak inter-triangle exchange 
(yellow bonds). At low temperatures each triangle behaves as 
a spin 1/2 entity (Sa = S' A = 1/2). Right: Energy diagram 
for one scenario of inter-triangle exchange where a term of the 
form (A/2)(SaxS' Az — SazS' Ax ) admixes several of the states 
\S,Ms;Sa = 1/2, S' A = 1/2 >, shown as blue lines. 



(shown in blue) are equal (J a ~ 65 K in units of ks), 
and an order of magnitude larger than the third (shown 
in red, J c w 7 K). Additionally, from nuclear magnetic 
resonance (NMR) studies and chemical structure anal- 
ysis it has been argued that there exists a very weak 
inter-triangle exchange interaction (yellow bonds, ap- 
proximately 0.3 K). In the absence of inter-triangle ex- 
change, for B — the ground state of each triangle 
consists of a 2-fold degenerate doublet with total spin 
S = 1/2, consistent with Kramers' theorem. The excited 
levels are a second degenerate doublet with S — 1/2 and 
excitation energy (J a — J c ) w 58 K, and a 4-fold de- 
generate level with S — 3/2 and excitation energy (also 
measured from the ground state) 3J a /2 w 97 K. In the 
experiments described below we consider temperatures in 
the range 1.5-5 K and B < 25 Tesla, well below the field 
value (w 74 Tesla) when the S = 3/2, M s = -3/2 level 
crosses the ground state S = 1/2, Ms = — 1/2 level. As 
such, it suffices to consider only the ground state doublet 
of each triangular unit. A weak residual inter-triangle 
anisotropic exchange will lift the 4-fold degeneracy for 
B = of each molecule and give rise (in general) to four 
distinct energy levels (see below). As remarked above, 
the occurrence of these splittings can be manifest when 
the molecules are subject to pulsed magnetic fields, giv- 
ing rise to a sudden reversal in magnetization when the 
field crosses B — 0, as a result of LZS transitions be- 
tween the split levels. Apart from the vicinity of B = 
the magnetic properties at low T of a {Ve} sample may 
be accurately described in terms of an ensemble of inde- 
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FIG. 2: Measured magnetization vs magnetic field for T = 1.7 
K (a) and T — 4.2 K (b), for the three half-cycle sweeps shown 
in the inset of Fig. 2(a). The solid red lines are obtained using 
Eqs. (1) and (2). The time dependence of the relaxation rate 
t~ according to Eq. (2) for T = 4.2 K is given in the inset 
of Fig. 2(b). The dotted line indicates the residual constant 
Ro in Eq. (2). The sudden drop of M to zero at B — is 
explained in the text in the context of LZS transitions. 

In Figs. 2 and 3 we present our experimental and 
theoretical results for the magnetization versus applied 



3 



magnetic field for two different temperatures (1.7 K and 
4.2 K) and for half-cycle and full-cycle sweeps shown in 
the insets of Fig. 2(a) and Figs. 3(a) and 3(b). The 
two striking features of the M vs. B data are hysteresis 
loops and the appearance of magnetization steps (in Fig. 
2) and near-reversals (Fig. 3) in the immediate vicinity 
of B = 0. The hysteresis loops (all data except in the 
immediate vicinity of B = 0) are reproduced (solid lines 
in Figs. 2 and 3) by numerical solution of the following 
generalization of the familiar Bloch equation^" 
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Here p, B is the Bohr magneton, p denotes the mass den- 
sity, v the sound velocity, V s i the characteristic modu- 
lation of the spin energy under long-wavelength acous- 
tic deformation, and M eq (T,B(t)) is the standard two- 
level equilibrium magnetization for an instantaneous field 
B(t) and for temperature T, i.e., M eq (T,B(t))/M max = 
t&nh[ gf i B B(t)/(2k B T)}, where M max = 2(Na£U B /2). 
We have derived Eq. (1) from first principles 10] upon 
making the assumption that for these temperatures the 
phonons are in thermal equilibrium with the cryostat at 
all experimental times. The first term of Eq. (2) is the 
low-temperature relaxation rate of the spins due to direct 
one-phonon processes, where spin flips are triggered by 
an acoustic phonon mode meeting the resonance condi- 
tion for the instantaneous energy separation of the two- 
level spin system. This term is a generalization of the 
standard expression for the relaxation rate due to one- 
phonon processes in a static external field 4] . Both the 
B(t) 3 factor, proportional to the phonon energy density, 
and the statistical mechanical factor depend on the in- 
stantaneous resonance frequency, proportional to B(t). 
The numerical value of V s i depends on the specific de- 
tails of the spin-phonon coupling (see, for example, dis- 
cussion for paramagnetic spins in Ref. Q), which at 
present is unclear. The quantity Rq in Eq. (2) rep- 
resents additional relaxation processes present and it is 
taken as a fitting parameter. Using the measured value 
of p — 1.93 g/cm 3 and estimating v — 3000 m/s, we 
obtain excellent agreement with our data for the choices 
R = 0.2 ms" 1 for T = 1.7 K and R = 0.5 ms" 1 for 
T = 4.2 K and V s i/k B = 0.35 K. Despite the small- 
ness of Rq it is important to retain this term in order to 
achieve a good fit to the experimental data in the low- 
field regime (below 4 Tesla for 1.7 K and below 7 Tesla for 
T = 4.2 K); for higher fields the dominant contribution 
to 1/r comes from the one-phonon term. We empha- 
size that the solution of Eq. (1) is extremely sensitive 
to the explicit functional form of the first term of Eq. 
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FIG. 3: Measured magnetization vs magnetic field for T = 1.7 
K (a) and T — 4.2 K (b), for the full-cycle sweep shown in 
the inset of Figs. 3(a) and 3(b) (dotted lines). The solid 
lines are obtained using Eqs. (1) and (2). The insets also 
show the measured magnetization (circles) vs time. The LZS 
transitions occur in the interval between the points A and B. 



(2)): Adopting a different choice of functional form one 
cannot achieve quantitative agreement with the observed 
hysteresis loops, for the whole field range and for differ- 
ent choices of field sweeps. Achieving an excellent fit to 
our measured data for a variety of choices of B(t) thus 
affirms the basic correctness of these two equations. 

We now discuss the magnetization steps observed for 
B « 0, and in particular the interval between points A 
and B in Figs. 3(a) and 3(b) (the other steps seen in 
Figs. 2 and 3 have the same physical origin and will 
not be discussed separately). In this interval the exter- 
nal field varies approximately linearly with time, with 
sweep rates of order 1 Tesla/ms. We note that at point 
A, MA/M max is somewhat less than unity due to ther- 
mal relaxation before entering the fast-reversal regime. 
Equivalently, at point A we are dealing with a statisti- 
cal mixture of spin-up and spin-down states. The most 
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striking feature though is that the magnetization Mb, 
at point B, nearly equals — Ma- We find that the time- 
widths of the near-reversals is shorter the faster the sweep 
rate, and is in the range 0.5-0.8 ms. We propose that 
adiabatic LZS transitions are responsible for the magne- 
tization steps observed in our system. The characteristic 
energy gap A of the LZS 2-lcvcl 13J model is related to 
the time- width St^zs of the magnetization step and the 
field sweep rate r by the relation StLzs = 2A/(<^sr) 
0,13). Thus, the measured time- widths of the steps give, 
as a first estimate for the zero-field energy gap, A » 0.4 
K. Using the above estimate for A, we are indeed in the 
regime of adiabatic LZS transitions, since the transition 
probability P LZS = l-exp(-Tr(A 2 /(2hgn B r)) « 1 0,0, 
thus implying that Mb = —Ma- The observed deviation 
of M/M max from exact reversal is about 15% for T = 1.7 
K. This discrepancy may be due to the role of the heat 
bath, i.e., the problem of dissipative LZS transitions (see, 
for example, |ll| and references therein); however, a sys- 
tematic investigation of this issue is in progress. More 
generally, it should be noted that, according to the above 
formula for Stj^zs, it is the high sweep rates used in our 
experiment that ensure that the adiabatic LZS transi- 
tions take place over such a short time interval as to be 
clearly distinct from the hysteresis loops of the thermal 
relaxation regime. 

The origin of LZS transitions in {Ye} remains to be 
discussed. The relatively large estimated value (0.4 K) 
of the zero-field energy gap A for these molecules suggests 
that its origin cannot be due to dipolar or hyperfine fields. 
In addition, as explained above, the lowest energy levels 
of each independent triangle are doubly degenerate for 
B = (and not four-fold degenerate as in {VisAsejQ) 
consistent with Kramers' theorem. Hence, for B = the 
ground state of such a molecule would consist of four de- 
generate states, namely the three symmetric states of the 
5=1 triplet and the antisymmetric 5 = singlet state. 
Inter-triangle exchange coupling could lift this degener- 
acy and give rise to avoided level crossings. However, 
since isotropic inter-triangle exchange cannot account for 
admixing states of total spin, we suggest that the avoided 
level crossings are due to the anisotropic (symmetric or 
antisymmetric) portion. One scenario is given in Fig. 1. 
The behavior of the dynamical magnetization thus in- 
volves LZS transitions between (at most) four levels. A 
detailed theoretical treatment of these transitions will be 
given elsewhere. 

In summary, time-resolved magnetization measure- 
ments using sweep rates of order 1 Tesla/ms show hys- 
teresis loops and magnetization steps for B w in the 
magnetic molecule {Vg}. The two effects are clearly dis- 
tinct because of the relatively high sweep rates used in 
our experiment. In the absence of both an anisotropy 
energy barrier and the phonon bottleneck effect, the 
hysteresis effects exhibited by this molecule occur be- 
cause the spin relaxation times are of the experimental 



time scale. Using a generalization of the Bloch equa- 
tion we were able to reproduce our experimental data 
for T = 1.7 K and T = 4.2 K for a large variety of 
field sweeps, and thus identify direct one-phonon reso- 
nant transitions among the Zeeman-split doublet of each 
triangle as the dominant mechanism underlying the hys- 
teresis behavior. The main assumption of our model, 
namely that the phonons are in equilibrium with the 
cryostat, should break down for temperatures below 1 K 
due to the phonon bottleneck effect. In fact, our prelim- 
inary data at T = 0.6 K indicate that the relaxation rate 
1/r deviates from Eq. (2). A systematic investigation 
of this issue is in progress. The steps of the magneti- 
zation for B » are attributed to adiabatic LZS tran- 
sitions between lowest magnetic energy levels impacted 
by the existence of anisotropic inter-triangle exchange 
interaction of order 0.4 K. This estimate is consistent 
with that previously suggested by NMR dataQ. A more 
precise value of A could possibly be determined by spe- 
cific heat measurements [1 2j . or by Electron Paramagnetic 
Resonance (EPR) techniques. The small departures from 
complete magnetization reversal suggests that one cannot 
entirely neglect the role of the heat bath. More gener- 
ally, exploring nanomagnets with pulsed magnetic fields 
can reveal a variety of fascinating dynamical phenom- 
ena and provide microscopic information that otherwise 
is not readily accessible. 
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